Detection of non-Gaussian Fluctuations in a Quantum Point Contact 
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An experimental study of current fluctuations through a tunable transmission barrier, a quantum 
point contact, are reported. We measure the probability distribution function of transmitted charge 
with precision sufficient to extract the first three cumulants. To obtain the intrinsic quantities, 
corresponding to volt age- biased barrier, we employ a procedure that accounts for the response of 
the external circuit and the amplifier. The third cumulant, obtained with a high precision, is found 
to agree with the prediction for the statistics of transport in the non-Poissonian regime. 



Recently, measurements of current fluctuations arising 
from charge discreteness have become an invaluable tool 
in mesoscopic physics, the most noticeable achievement 
being the shot noise measurement of quasi-particle charge 
in the fractional quantum Hall effect [H, 13] • Typically, 
mesoscopic shot noise experiments report zero-frequency 
noise power, but this quantity contains only partial in- 
formation about the statistics of the transmitted charge. 
The counting statistics (CS) [3( is entirely characterized 
by the set of cumulants (irreducible moments) ((<?")) of 
the charge q(r) transmitted through a voltage-biased sys- 
tem during a sampling time r. In the long time limit 
they are proportional to the time, ((<?")) = ((J n )) t; this 
expression defines the current cumulants ((J™)). For ex- 
ample, Gaussian noise is fully determined by the first 
two current cumulants, the average current I = ((J)) and 
the noise power S = ((J 2 )}. The simplest measure of the 
non-Gaussianity, the third current cumulant S — ((J 3 )) 
which reflects the skewness of the current distribution, 
is the central focus of our paper. It is linear and uni- 
versal at low bias voltage, and therefore may be used as 
a tool for investigation of strongly correlated systems, 
where large bias cannot be applied without substantially 
affecting their properties. 

However, during a typical sampling time a large num- 
ber of electrons passes through the system. This fact, by 
virtue of the central limit theorem, makes it difficult to 
observe non-Gaussian effects in electron transport [4|], un- 
less electrons are counted one by one as, e. g., in Coulomb 
blockaded quantum dots [y, |7|, l8( . To date, ((J 3 )) has been 
measured only in low transmission tunneling junctions ei- 
ther by measuring voltage on a load resistor 0, [nj, or 
with the help of the on-chip Josephson threshold detector 

n 

In a typical experiment, due to the electrons' charge, 
the voltage across the sample is not constant, so the 
measured statistics are not trivially related to the CS 
of the voltage-biased system. Indeed, the original ex- 
criment on the third cumulant of a tunneling current 
exhibited totally unexpected results, which were ex- 
plained EH, 0, El by the back-action of the mea- 
surement apparatus on the sample. That is, in addition 



-1.3 -1.2 -1.1 -1 

Gate voltage [V] 

FIG. 1: The conductance step of go = 2e 2 jh versus gate volt- 
age V g used to control the transmission F of QPC conduction 
channels is shown for one of our samples. Lower inset: The 
Fano factor F2 = S/el versus the QPC conductance G/ go for 
the same sample. Solid line is the theoretical prediction, Eq. 
0, for the large bias, eV > k B T (T = 1.8 K, B = 2.2 T). 
Upper inset: The Fano factor F3 = ({J 3 ))/e 2 I versus V at low 
bias, eV « k B T (T « 5K), for QPC1 (♦), and QPC2 (■). 
Theoretical prediction F3 = (1 — V) is shown by the solid line. 



to ((J 3 )}, the experimentally measured potential fluctu- 
ations also contain contributions, dubbed "cascade cor- 
rections" in Ref. [l3| , from the voltage dependent current 
/ and noise power S. Surprisingly, even if the load re- 
sistance is small, the cascade corrections may be of the 
same order as ((J 3 ))- 

We report here the first measurement of the third cu- 
mulant of the non-Poissonian current partitioned by a 
variable transmission barrier. To obtain the third cumu- 
lant, we develop a procedure that allows us to separate 
the ((J 3 )) contribution and cascade corrections in a re- 
liable fashion. In particular, the frequency dependence 
of the amplifier gain is found to affect these contribu- 
tions differently. The resulting cumulant ((J 3 )} agrees 
accurately with the predictions of [3| for all transmis- 
sions, < r < 1, without fitting parameters. As a 
variable transmission barrier we use a quantum point 
contact (QPC) - a small quasi-lD constriction formed 
in a 2D electron gas by negative voltage V g applied to 
split gates [id] . By varying V g one can gradually in- 
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FIG. 2: Experimental setup schematic. The current source 
drives constant average current I through the QPC sample, 
while the capacitor CV fixes the voltage across the sample. 
The QPC resistance R s can be tuned by the gate voltage V g . 
Ri is the load resistance and C s t is the stray capacitance of 
the wires at the amplifier input (the point N). Current fluctu- 
ations generate voltage fluctuations at the point N, which are 
amplified and digitized by 12 bit analog-to-digital convertor 
(A/D). We introduce current through the capacitor C c and 
measure the response with a network analyzer to calibrate the 
setup. 



crease the width of the constriction. Fig. [T] shows a typ- 
ical dependence of the QPC conductance on V g . The 
go = 2e 2 /h step in conductance vs. V g , and the ac- 
companying noise minimum (bottom inset), defines the 
lowest spin-degenerate QPC channel used in this work. 
We define the average QPC transmission as: 



r(/) = 



9oV s {iy 



(1) 



where / is the average current and K is the voltage across 
the sample. The observed noise agrees well with the the- 
oretical expectations [l7| and has only thermal and shot 
noise contributions. No material contribution, e.g. 1/f 
noise, is detected. The measurements are done at el- 
evated temperature T « 5 K in order to reduce nonlin- 
earity of the current- voltage characteristic. Nevertheless, 
for QPC conductance below 2e 2 /h, the transport is dom- 
inated by a single channel, as confirmed by the agreement 
of the measured noise with the single-channel expecta- 
tions (see Fig.[3|). 

It is instructive to discuss the issues encountered in the 
investigation of the third cumulant prior to presenting 
our experimental results. We focus on the experimental 
setup shown in Fig. [2l however our conclusions are quite 
general. Current fluctuations J generated by the QPC 
sample and by thermal noise in the load resistor Ri give 
rise to the voltage V/vM = Z(u)J{uo) at the amplifier 
input (point N in Fig. [2]), where the impedance Z(lu) = 
R\\/(l — iu}TRc) is determined by i? 1 1 = R s Ri/(R s +Ri) « 
7KOhm, and by the time constant trc — R\\C\\ (here 
C|| = C st + C c « 3pf) which sets the high frequency 
cut-off of the circuit 1/(2tttrc) ~ 7 MHz. The amplified 
voltage is: 



KM = K(w)Vn = A(u)J(u), 



(2) 



where K is the amplifier gain and Vjsr is the input voltage. 

In the long time limit (which is justified since the 
high frequency cut-off of the system is much smaller then 
max(eV s ,kBT)/h) all current cumulants are frequency- 
independent. The amplified voltage fluctuations ((K 2 )) = 
([SV a (t)} 2 ) bear a simple relation to the noise power 



((V 2 )) = B 2 S, B 2 



1 

2^ 



A{w)A{-w)dw. (3) 



In contrast, the third cumulant ((K 3 )) 
be decomposed as 



[SV a (t)} 3 }, may 



<(K 3 » = B 3 ((J 3 )) + B ca J cn + B al J nl 



(4) 



That is, apart the third current cumulant ((J 3 )), it also 
contains "environmental" cascade correction J cri origi- 
nating from the back-action of the voltage fluctuation 
across the load on the current fluctuations in the sample 
at later times [l2, 13, 14, 15 1, and the correction J n \ due 
to nonlinearity of the sample: 



J en = 3R {l SdS/dV N 
J nl = 3RlS 2 d 2 I/d 2 V N , 



(5) 
(0) 



where S — S s + Si is the noise power generated by the 
sample and the load. The coefficients B$, B cn , B n \, de- 
rived below, depend on the circuit only and are of the 
same order. The thermal noise of the load and its re- 
sistance are current-independent, and thus do not con- 
tribute to the derivative. Note that since R\\S — ► 2ksT 
at — ► 0, reducing load resistance does not eliminate 
the corrections (0 and ((6|). 

Our measurements are performed using a setup (see 
Fig. [2]) similar to the one discussed in detail in [lflj. We 
improved it by increasing the total gain of the amplifi- 
cation chain to utilize the full 12 bit resolution of the 
A/D converter, thus eliminating the necessity to com- 
pensate for its nonlinearity We also managed to reduce 
the nonlinearity of the cryogenic amplifier by increasing 
the current through the transistors. The amplified signal 
is computer analyzed to construct the probability dis- 
tribution function (PDF) of the amplified voltage (|2]). 
We calibrate the setup by introducing a known ac sig- 
nal, spanning the entire frequency band of the amplifier, 
through a small capacitor C c « 2.4 pF. This gives us the 
complex-valued A(u>) [see Eq. @] independently at each 
value I and V g , which is subsequently used to compute 
the coefficients in Eqs. <j3j) and |4|. 

To compensate for inaccuracy in the calibration and 
for the drift of amplifier parameters we multiply the mea- 
sured A(u>) by a numerical factor, which scales the mea- 
sured noise S to the theoretical prediction [17[ 
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FIG. 3: The current cumulant ((J 3 )} at Y ~ 0.3, derived from 
the experimental results for ((K 2 » and ((V a 3 )} using Eq. (|4]|6)l 
for the sample QPC2 measured with amplifier b. Solid line is 
the theoretical prediction for ((J 3 )), Eq. ((HJ ; the dashed line 
is the low-bias limit, Eq. ((9j), with current-dependent trans- 
mission r (lower inset). Upper inset: Measured value of the 
noise power S versus current; solid line is the Eq. {7J|- 

This scale, determined individually for each value of V g , 
is found to deviate from unity by less then 10%. This 
is illustrated in the upper inset of Fig. [3j which shows S 
fitted by Eq. J7|) at V g corresponding to the transmission 
r « 0.3 (see the lower inset of Fig. [3]). We use the scaled 
A(oj) to find the coefficients B3, B en , B n \, and then to 
obtain ((J 3 )) from Eqs. |4} and (©. The resulting ((J 3 )) 
is shown in the main panel of Fig. [3l 

As seen in Fig. [31 the measured third current cumulant 
((J 3 )) shows very good agreement with the prediction for 
noninteracting fermions [H, Q : 



<(J 3 ))=e 2 /(l-r) 



6r smh(LQ - U 
~U cosh([/) - 1 



1 - 2Y 



(8) 



Although the procedure that leads to ((J 3 )) involves sub- 
traction of several terms of comparable magnitude, we 
stress that it does not rely on any fitting parameter other 
than the aforementioned scaling factor. Since we observe 
no systematic deviation from the prediction ([SJ, we be- 
lieve that the main sources of error in our experiment are 
statistical fluctuations (indicated by error bars in Fig.[3j, 
as well as our lack of knowledge of the precise energy de- 
pendence of r. 

The nearly linear behavior of ((J 3 )) at I < 10 nA cor- 
responds to the low-bias limit of Eq. ([5]) : 



((J 3 )> = e 2 /(i-r), \u\«i. 



(9) 



The dashed line in Fig. [3J shows the prediction © with 
the measured current-dependent transmission r(/). No- 
tably, the full expression ([5]) provides a much better fit 
to our data than the low-bias limit ([9]). In the data taken 
at larger bias (see Fig.[4j, for which the quantity U could 
be as large as 4, the agreement with the expression (|8]) 
remains very good. We also note that the large bias limit 



of Eq. 



((J 3 )) = e 2 I{\ - T)(l - 2Y), exhibits a sign 
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FIG. 

Eqs. (Hi): <(V a 3 »/B 3 (■), JenB en /B 3 (*), and J nl B nl /B 3 
(♦). Panels (a) and (b) show the results obtained with two 
different amplifiers, a and b, on the sample QPC1 for similar 
transmissions V — 0.6 and 0.7. The coefficients Bz, B cn , and 
_B n i are calculated using measured response function A(ui). 
Solid line is the prediction p). 

to observe this effect, U > 10, is not accessible in our 
experiment because of the nonlinearity in J(Vg). 

In this work we present the data obtained on two dif- 
ferent samples, QPC1 and QPC2. The sample QPC1 is 
measured with two amplifiers, a and b, which have dif- 



ferent low frequency cut-offs, i/< ~ 300 KHz and 
4 MHz, and upper cut-off frequency z/> 3> 1/Q-ktrc) ~ 
7 MHz. For comparison, we show the results for two sim- 
ilar QPC transmissions in Fig.|4j Although the contri- 
butions to ((J 3 )) are quite different, the resulting current 
cumulants are close and agree well with the theory. This 
underscores the necessity of accounting for the frequency 
dependent amplifier gain A(uj) in data processing. The 
low bias results for two different samples are summarized 
in the upper inset of Fig.[TJ the slope of the ({J 3 )) vs. cur- 
rent, the Fano factor F$ = ((J 3 ))/e 2 / is proportional to 
(1 — T), confirming universality of the relation ((9]). 

In the rest of the paper we derive B en and B n \, and 
express them through A(lj) and the input circuit param- 
eters. Fluctuations are treated in the time domain us- 
ing the Fourier transforms A(t), Z(t), and K(t), which 
vanish for t < due to causality. The fluctuation 
V(t) = Vn — IRi of the voltage VnH) is generated by 
the current noise J according to 



CuV 



-Rn X V ■ 



d 2 I V 
dV 2 2 



J. 



(10) 



Solving linearized equation (|10p yields an exponential re- 
laxation of fluctuations: 



V(t) 



db\Z{t - ti)J(ti), 



(11) 



change at Y = 1/2. However, the bias regime needed 



where Z{t) = Cm 1 e~ t l TRC . The amplified voltage V a is 
related to V(t) as: 

V a (t)= [dt'K(t-t')V(t')= fdtiA(t-ti)J(ti). 



4 



We next ensemble-average the value [V a (t) — (V a )] 3 . 
The intrinsic contribution arises from (( Jfa) Jfa) Jfa))), 
where Jfa can be treated as ^-correlated in time. This 
leads to the first term in Eq. (H|) with B 3 given by [lfj] : 



(2^ 



+ 00 



A(u)A(u')A(-w - J)dujdJ (12) 



The voltage across the sample fluctuates in time, be- 
ing a function of current fluctuations at preceding 
times, as described by (jlip . As a result, the average 
(Jfa) J fa) J fa)) contains the environmental contribu- 
tion due to the dependence of the correlator 



((Jfa)Jfa))) =Sfa -h)S(V) 



(13) 



on the voltage V, which in turn depends on the current 
J at an earlier time t'. Linear in V expansion of Eq. (JTJ 
gives the environmental contribution to Eq. ((4]) : 



3 / dhA 2 (-ti) 



dS 
dV^ 



dt'Zfa - t') 



x{{J(t') / dt 3 A(-t 3 )Jfa))) = 3R l{ S 



dS 



-B c 



B r 



R» 



dhA^i-h) / dt 3 Z(t 



dV N 
t 3 )A(-t 3 ), (14) 



where the factor 3 accounts for the three possibilities to 
choose a later time. 

Finally, the nonlinear contribution to Eq. ([4]) comes 
from the V 2 term in Eq. (fT0|) with the coefficient 



B nl = Rj 2 J dtA(-t) (^J dtiA(-ti)Z(t-ti 



(15) 



For the frequency independent amplification K, we find: 

B 3 = 2B ca = 4B nl = K 3 t rc /3C 3 . (16) 

In this case, the first two terms in Eq. ((4]) cancel at small I 
and Ri. Our amplifier a, having v < ~ 300 KHz, operates 
close to this regime (see Fig. SJ). 

It is instructive to note that the ratio of the intrin- 
sic and environmental coefficients in Eq. (fT6|) is twice as 
large as that of Refs. [II El, H, [II]. This difference 



can be traced to different assumptions about frequency 
dependence of the system gain. Refs. [H Q 111 fo- 
cus on the limit of the high frequency cutoff set by the 
amplifier, while in Eq. (| 16[) . as well as in our experiment, 
it is determined by Z(u) with the roll-off set by 1/trc- 
Therefore, our results correspond to the equal time cor- 
relator ((V 3 )) = {[8V{t)] 3 ). 



In summary, we have measured the third cumulant of 
shot noise in variable transmission QPCs in the essen- 
tially non-Poissonian regime. In order to extract the "in- 
trinsic" third current cumulant we developed a technique 
which allows reliable subtraction of the environmental 
and nonlinear circuit-dependent contributions. Good 
agrement between the experimental results and the ex- 
pectations, Ref. [H , opens a venue for using high order 
cumulants as an experimental tool in mesoscopic physics. 
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